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ABSTRACT 


SA This research investigates seismic signal processing 
techniques for battlefield target classification and acgui- 
sition. Multiple target classification is performed by 
Giscrete time domain matched filtaring. Multiple ‘target 
directions are determined using the responses of the matched 
filters andleast sean squares polynomial curve fitting. 
The least mean squares polynomial carve fitting procedure is 
also used for direction finding for recoil/blast scurces, 
using the unfiltered seismic signals. 
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I. LN RODUCTIIN 


Timely end accurate combat intelligence is a integral 
part of the aodern battlefield. A primary goal of combat 
intelligence is target acquisition. Hostil2 targets can be 
acquired by either passive cr active aeans. An example of 
passive target acquisition is visual target ilentification. 
Radar, on che other hand, is an active target acquisition 
device. An effective combat iate2zlligence system will 
include a mix of both active and passive target acquisition 
methods. 

Rapidly aivancing technology in the fields of electronic 
counter-measures and radiation-seeking weapons has enhance? 
interes< in passive target acquisition mecieds. To be cos* 
effective, as an additional target acquisition system, a2 
passive system must be able to proviie swift and accurate 
target identification, location and tracking information on 
hostile targets. A variety of seismic sensor systems have 
been used in this roll with varving degrees of sucess. 

Naval Ocean System Command (NJSC) in San Diego has 
developed a system based on a circular ring c£ sensors with 
data collection managed by an arcay-proecessor/minicomputer 


sys<em (Ref. 1]. The observation of enemy movemerts and 
activity beyond the Porward Edge of tha Battle Area (FEBA) 
is the design objective of this system. This thesis uses 


data collected during a test of this system at the Marine 
Corps Air-Ground Combat Center at Twenty-nine Palms, 
Califiornia. Investigated are varidsus methods of processing 
the seismic data collected. The objective of this research 
is to try to provide viable methods of satisfying system 
design objectives through signal processing techniques. 
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Tha following chapter addresses the design objectives 
and requirements for such a seismic Sensor System. 
Additionally, capabilities and d2ficiencies of current 
systems and research are detaiied. In order to intelli- 
gently address solutions to these r3quirements and deficien- 
cies, an understanding of seismic theory ani sensors is 
needed. Genaral seismic theory is presented with emphasis 
on the constraining parameters for the use of the earth's 
surface as a medium for gathering saismic intelligence. 
Also highlighted are the similaritias of the earcth's surface 
to electro/optical phenomena and the resulting simplifying 
assumpticns. 

Tha sequence of design solutions investigated followed 
from the analogies and simplifying assuaptions addressed ir 
the study of seismic senscr theory. The problem cf targe+ 
identification or classification is aoproached using digital 
matched filtering cf the time domain amplitude data. 
FPraguency domain matched filtering was not considered, based 
upon the conclusion by NOSC that there appeared to be no 
consistant spectral lines for any of the possible target 
types, except for artillery [Ref. 1]. Matched filtering was 
used to identify single and sultipla target classes occuring 
during a Sample period. The chaptar detailing the matched 
filter procedures and implementation also includes 3 
description of support and walidiation software used in the 
analysis of the seisagic data. 

The validation software is priaarily used to check the 
accuracy of the direction finding routines. These direction 
finding routines are the tine domaia phase difference proce- 
dura (TDPD) and a least mean squares polynomial (LMSP) curve 
fitting procedure. The combining of the matched filtering 
procedure and these direction routines allowed for auitiple 
target direction finding. The theory and derivations of the 
two direction finding algorithus are presented in the 
aultiple target direction chapter. 
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Application of these algorithms is performed first on 
Simulated targets for validation of the procadures. The 
experimental data is then analyzed. A user's manual, which 
includes procedures for tape and mass storage operations, is 
provided as an appendix. This appendix describes how te set 
up and use the software system. 
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TI. THEORY OF SEISSES SENSORS 


Elastic waves result from the stressing o£ ar elastic 
media. The alastic media for seismic theory is the earth. 
Seismic thecry is the study of tha earth as a wave propa- 
gating media. Elastic waves propagate away from the source 
of seisnic stress, e.g. an explosion {[Ref. 2]. The energy, 
which propagates through the earth, travels via particle 
defermations. The elastic properties and densities of the 
earth media determine the velocities of these seismic waves. 
(Ref. 3]. Seismic wave sources of interest may be inpulsive 
or continucus. Impulsive or short duratica sources ar? 
artillery recoil or shell blasts. The time-limited nature 
of this type of seismic signal produces a broad range of 
frequencies. Continuous wave signals may be produced by 
tanks, trucks and lew flying aircraft. These continuous 
wave signals may be described by narrow band frequency char- 
acteristics. The spectral power of a seismic source is a 
function of several parameters. A non-inclusive list of 
these parametars includes: 

1. The vehicle's velocity and a2ass 

2. The size of the explosive charg2 associated with the 
artillery or sheil blast 

3. The degree of coupling into the earth's surface 

4. The geological structure over the wave's path 

Information about the seismic source is contained in the 
waves which it generates. Por axaaple, in an array of 
seisnic sensors, iirectional information is contained in the 
relative received signal phases. In otherwords, the rela- 
tive phase differences between the signals received by the 
individual array elements can be used to compute the direc- 
tion to the saismic source [Ref. 4%]. These relativs phase 


differences represent the time delays of the waves as thay 
pass the array's elements. The response of the arzray's 
elaments is proportional to the aaplitude and velocity of 
the earth's motion, relative to tha geophone's sensing axis 
fer the waves. (Ref. 3] 

Assumptions about the propagation of seismic waves must 
be made to assist with their analysis. The earzh is assumed 
to be made up of horizontal, homogeneous, and isotropic 
layers of material. These layers ara assumed to be discon- 
tinuous in their alastic properties at their borders. This 
variance in the elasticity between the layers leads to an 
optical analogy for the wave propagation. Prcpagation paths 
may now be viawed as being direct, refracted, or reflected 
versions of the source's seismic waves. (Ref. 3] 

There are four Lasic types of seismic waves. These 
types are compressional, shear, Love, and Rayleigh waves. 
Compressional waves are generated by impulsive sources such 
as shell blasts. Particle motion is along the direction of 
travel. Shear waves are charactacized by particle action 
orthogonal to the direction of motion. The Love wave is a 
surface wave which may occur as a result of the layering of 
the earth's surface. This layering effect acts as a wave 
guide for this type of wave. Particle motion is orthogonal 
to the direction of wave propagation. The Rayleigh wave is 
generally the strongest of the seisaic waves. The Rayleigh 
wave travels along the free surface of the earth. Its 
particle motion direction is always in the vertical direc- 
tion. It is the strongest wave generated by a compressional 
source. Its aaplitude attenuates at a rate only inversly 
proportional to the square root of the distance. [(Ref. 3] 

The waves of primary interest for a seismic system are 
the Rayleigh and Love waves. This is due to the long-range 
propagation of these waves. Since these two wave types are 
orthogonal to each other, they must be sensed by different 


tN 


2 TRS 


geophones. Rayleigh waves may be sensed by vertical 
geophones and the Love waves by horizontal geophones. Since 
the Rayleigh wave is normally the strongest, and in order to 
reduce computational complexity, only verticle sensor iat3 
is used. 

Rayleigh waves experience absorbtion losses, particu- 
larly at higher frequencies. This phenomenon occurs because 
of the lowpass filter effect of ths earth. This filtering 
effect is further compounded due to the fact that the cutoff 
frequency of the earth diminiskes with range. Purt her 
complications arise due to the dependence of wave velocity 
en frequency. The result being that che wave “rain may 
change with distance, reducing the correlation of the wave 
shape between its source and distant points. Other sources 
of error occur because of the weathering of the surface 
layer, irregularities in the sub-surface composition, vari- 
ances in the earth's layers and surface geometry. (Ref. 3] 
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TIX. DEVELOPHENTAL BEQUI REM ENT 


A. PASSIVE TARGET ACQUISITION AND SURVEILLASCE 


A battlefield commander possess2s a definate raquitement 
for real time combat intalligence. A significant tactical 
advantage is held by the commander who is able to intecrate 
his availiable combat intelligence sources with his 
supporting arms, i.e., ‘target acyuisition and engagement. 
It follows that to be part of ths target acquisirion 
process, any ceal tine, seismic sansing syst2m must provide 
swift and accurate information on datected 2nemy targets. 
Tha specific requirements for such a system are the ability 
to detect, identify, and locata these targets [Ref. 5}. 
Additionally, the target's rate and direction of movement 
should be provided or made easily discernible. 

Any seismic sensing system must be designed around the 
target acquisistion cycle. The target acquisition cycle, as 
given by Dublin (Ref. 5], is as follows: 

1. Search Time 

2. Target Sensing 

3. Information Processing 

4. Display of Target Information 

5. Analysis of Target Information 

6. Time required to make a Decision 

7. Time Required for Supporting Weapons to Resvcend 

For a seisaic system, a prioritized list of possible 
targets are as follows: 

1. Artillery 

2. Helicopters and Aircraft 

3. Tracked and Wheeled Vehicles 


4. Perscnoel 
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As may be expected, the relative amplitudes of these seisnic 
targets vary widely. A seismic targeting system is thezre- 
fore constrained as to the targets it can or can snot be 
expected to effectively engage. 

The variance in the relativa azplitudes of saissic 
targets suggests a range of specifications for detection of 
thase targets. As summarized by Dublin [Ref. 5], possible 
detection radii may be as shown ir Table I. Radii are given 
for both short and leng targeting syst2us. 


poo 
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| TABLE I | 

: { 

| Target Detection Radii | 

! Target Short Bange Syst2a Long Raage Systen 

| Personnel 100M None 
Vehicles : 1KS 10-20KM 

{ Low Plying Aircraft 1K 1O-29KM | 

15-20KM | 

4 


| Hostile Weapons 1K 


The timeliness requirament is ancillary to the radii of 
datection specifications. Timelinass, as used here, refers 


to the total time coagsgencing wher thea seismic sensor system 
first detects and processes the seisaic target data and 
ending with the dissemination cof tha targeting information 
to command slesents for disposition. This tiseliness 
requirement ranges between five to fifteen ainutes, 
depending upon the sobility of tha target (Ref. 5). 

The parameters having a direct effect on the tiseliness 
of a system are the probabilities  2f false alara and detec- 
tion. These parameters directly relate to a system's value. 
Increased probability of detection with reasonable false 
alarm performance, coabined with tha ability to disreguard 


friendly targets, are practical design objectives for any 
targeting system. The ability to incorporate such design 
features into a seissic sensor system will reduce both thea 
time wasted on invalid targets and the danger of undetected 
*argets. 

Once a valid target has been detected, ‘arget location 
information must be obtained. Stringent specifications for 
target locations allow for the system to support or enhance 
the effectiveness of; fire support systems, blind bombing, 
Harassing and Interdiction fire (H and I), and observerless 
artillery ergageuen+. 


B. CURRENT CAPABILITIES AND DEPICIENCIES 


To date, numerous successful aijgorithms have beén devel- 

oped to determine direction tc single targets. These algo- 
ithms include both tima and fraquancy domain methods. 

Target identification of long range targets via seisnic 
sensing has not yet met with equal success. 

The modera battlefield is seldom a single ‘*azget types 
environment. The complicated, real woerli problems of 
multiple target identification ard multiple target engage- 
men+ require soluticn before practical seismic sensor 
systegas can be integrated into the target acquisition 
process. 


ra 


IV. HATCHED PILTER CONCEPT AND DISCRETE ALGORITHS 


A. UATCHED FILTER THEORY 


As previously addressed, there exists a reguiremens for 
battlefield target identification/classification. The | 
recovery and classification of targat signals suggests a 
filtering requiremert., Previous works ard implementations 
have used frequency domain technigues [Ref. 4]. The Air 
Force's SKEET system and the U.S Army's Remote Battlefiel3 
Surveillance System (REMBASS) both have successfully impie- i 
mented a spectral power approach for classification of 
seismic data. These systems, however, re for shert range 
applications. Time domain approaches +o target identifica- 
tion have been for the most part laft unexplored. 

The discrate matched fil~ering technique is an attemp+ 
to classify targets by thair time jomain amplitude pattern 
i.2., *heir seismic amplitude signature. The heuristic 
basis for this method evolved from the observation cf visual 
differences batween the amplitude versus time signals for 
the various classes of targets. The matched filter, being 
tha optimuga filter for detecting known signals, was selected 
(Ref. 6]. 

The discrete matched filter is described by equation 
4.1, where h(t) is the impulse rasponse of a filter whose 
output signal to noise ratio at tims t, is maximized. The 
uniz step u(*) has been added to assure causality for the 
system. The gaatched filter is {Ref. 7]. 


h(t) = s(*, -t)a(t) (4.1) 
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The output signal is given by 
oo 
scut - frome: - v)dv (4.2) 
— OF 
The maximum signal to noise ratio is given by 
(sout2?/M2)maxoutput = E(t, )/N% (4.3) 


Whare Mis the noise level at the filter outout, N is the 
input noise level, and E(t,) is the energy in s(t) up to 
time t,. In equation 4.2, the replica of the oriqinal knewn 
Signal is raversed and translated in time to be convclved 
with the signal input to the filter, producing the optimun 
output signal to noise ratio. 

Fer discrete realization of the matched filter as inple- 
mented, equation 4.2 becomes equation 4.4, where h(k) is the 
reflected and translated known signal. 


N 
sout (4) = Pa ckis(3 - k) (4.4) 


Where N is the number of data points per sample period. 


B. SATCHED FILTER ALGORITHSA 


Samples of known signals are stored in a data file and 
Tead into a5120 array at the start of program execution. 
Five sample, or known filter signals, are recorded in this 
array. Each sample signal comprisas 1024 9£ the 5120 
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elements of the matched filter array. To perform he 
matched filtering, the 1024 seismic data samples (1nknown 
signal) are copied cver five tims on an array of 11264 
size. These copies cf the experinental data are separated 
by 1024 zeros on aither side. Additionally, «he leading and 
trailing elements of the axperimental data are set to zero 
te gliminate switching spikes oresant in the data. Farce 
forth, this 11264 element array will be referred to as the 
working array for brevity. This working array will contain 
tha results of tha matched filtaring. For M, given as tke 
number of input signal and filter signal arsay elements and 


ag 
i 
‘s) 


also tke number of zeros, the requirement set fert!} 
{Ref. 8} for nonoverlapping convolution of length L is 
satisfied. Equation 4.5 establishes the minimum length for 
nonoverlapping convolution of one matched filter segment. 


L = 2M = 1 (4.5) 


Figure 4.1 depicts the working array layout. Notice 
that the first filter signal, "Trackad Vehicle" is lcaded 
into the h(t) artay and is convoluted with the working array 
elements one through 2048. Once the leading data point of 
h(t) reaches the working array's element 2048, a new known 
signal is loaded into h(t), “Wheelad Vehicle", and the 
convolution is continued for working array elements 2048 to 
4096. This process continues until the last of the five 
h(t) filter signals has baen read in and the convolution has 
been performed on all copies of the data. Notice that 
wotking data array elements 10240 to 11264 are for array 
length over-run protection. 

Prior to sonvolution, each filter sample signal and the 
input signal are equalized to the same power level through 
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i | 
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I { 
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* * 
. f { 
| { 
; I - Indicates point at which new xnown sample sianal | 
| is loaded into h.t}. ( 
I 
{ { 
{ | 
Arcee eens 
on _ —— ee | 
Pigure 4.1 Working Array Configuration 
division by their respective root mean square values. 
Additionally, after the convolution, the entire working 
array is normalized with respect to itS maxinun ampoliztude 
element. Tha maximuan val 5 


u2 in ease target classificatios 
section is then compared with tha interactively selected 
matched filter thresheli, If the section's peak values is 
above the threskold value, that class of target is declared 
to be present. This mathod allows ‘for the simultaneous 
detecticn and classificaton of multiple targets. It should 
be noted «hat this equalization forces beth high and low 
amplitude signals to an equivalent average power level. 
This is felt to be justified in that signal amplification is 
mot the goal, rather target detection and classification is. 
In the case of a single target, the known signal that most 
closely matches the unknown signal is anticipated ‘to have 
the greatest aaplitude matched filtar spike. Figure 4.2 is 


a sample output. 
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C. SULTIPLE TARGET BATCHED FILTERING 


The upper Lizait of equation 4.4 aay be selects? “=c 5s 
fron one to 1024 when cailed by the multiple target irec- 
tion routine. This allows for the selecrti be 
program execution times. Target iitnctifica 
always made oy the full 1024 elenenc bu 
matched filter target idantification routine is used by the 
multiple target aisecticn findin coutine, data windows of 
less than 1024 are formed from 2 seynent 29 the 1024 
elements by axtracting tis seca uated arcuné the 
Baxinum signal value. Figure 4.3 illustrates bis 


procedure. 


INPUT 


TT ay OT en EE Ee OS A ce OO aD Ate OS eee a me 


9 10 1 {2 13 4 5 16 17 


TIME ~SECONOS~ | 


Figure 4.3 Windowing of Experimental Data 


| 


D. COMPUTATIONAL REQUIREMENTS 


The computational requirements for this proceduse are 
Gescribed by aquation 4.6, where NJ is the “otal number of 
operations zsequired by the matched filter routine. 
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a EN eee oe ted 


NO = (Ndf2) (Nc) * (Ndf - 1) (%4£) (Nc) (4.8) 


f is the number of data elenents ana also the nunper of 
{ filter elements. Ne is the total numbe2r of tar 
Tha term Ndf2Nc is 
required, while (Ndf - 1)NdENc is the tctai a 
a 


get classes. 


a 
the total nuaber sf wulrtislicaticns | 
u 


mper of 


When used for target identificacion, Ydf aquals 1024 and 
Ne equais five. Equaticn 4.6 givizs a2 total of 10,439,640 
operations for thes as. ahen the matched 
routine is initiat 
routine, the value 
to 1024. Figure 4.4 is a plot of aguation 4.6 
computational consequence for selection of ia 
Ndf. Note *hat equation 4.6 gust 32 aultipnled 
or sensors in the ring when computing zhe nu 
tions for the multiple direction routine. 
: Window sizes 9f 100 and 200 were found to provide accep- 
: table accuracy with greatly reduczi computation tines. As 
shown by figure 4.4, a window sizs of 200 requizes 399,000 
operaticns, while a wiandsw size of 109 requires only 99,500 


operations. The number of operations required was found, as 


— 


expected, +o be proportional to the exscution time of this 
routine. 


E. SUPPORT SOFTWARE 


1. Amplitude Analysis ( Timout ) 


A gtaphical output is provided by this amplitude 
routine which displays the relative amplitude versus time 
for a selected sensor. The initial target direction and 
target classifications found are also displayed. This 
routine allows for the interactive selection of the 
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Figure 4.4 Data Window Size vs Number of Operations 


amplitude response of any sensor as a sample target for 
later use in the matched ‘filter analysis. In this way 
Sample, signals can be catalogued and evaluated as filter 
Signals. The axes of tha graphical output adapts to the 
data's maximua aauplitude and to the time period involved, 
Figure 4.5 is an example of amplicuie analysis output. 


2. Preguancy Aag.ysis ( Freget ) 


The frequency routine, as in the amplitude analysis 
routine, allows for initial primary target direction, target 
classifications, and for the time pariod of the experimental 
data to be displayed. Normalizai spectral power versus 
frequency is graphically displayed as shown in fiqure 4.6. 
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3. Simujacion and Validation Software (SIMULT) 

In order to valiijate tha various algorithns and 
their implementation in software, a testing procedure was 
required. The specific algcerithas which «he saimulatior 
routine was Jasigned to validat2 are the initial asgle, 
phase differance and least mean square target direction 
routines. These routines, as previously described, us2 «he 
ralativs tima differences of the seismic signal's peak 
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amplitude response cr the peak aatcned filter response 
respectively. Validation of these routines is performnsd by 
allowing the creaticn of simulatai ‘targets with s2lected 


arrival angles. 

Simulated ‘¢argets can be created in the routine 
SIMULT. Up to four sine wave targets of selected frequency, 
amplitude, and direction can be input during each sanpis 
period. These simulated targets are added to the experi- 
mental seismic signal data for the sample pericd. 
Correspondingly, zere diraction filtar data nust have been 
written into the matched filter data file at these selecte4 
sample target frequencies. 

The directions for the simulated targets are created 
by introducing relative phase delays batweern the sine waves 
that are addad to each sensor's seismic data. This is 
implemented by introducing a zero phase to the sensor in the 
desired direction of the simulated target. The phases of 
the other sensors are increased proportionally by <‘«heir 
distance onthe circular array away from the zero ohase 
shift sansor. Pigure 4.7 illustrates the case of a zero 
degree simulated target direction. 

A summary of the target dJiractions found by the 
multiple target direction routine and the simulated targets 
entered is provided by the aultiple direction plotting 
routine. Figure 4.8 is a sample output. 
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Simulation of a Zero Degree Target 


Figure 4.7 
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VY. MOLTIPLE TARGET DIRECTION 


A. THEORY AND DESIGN CRITERION 


The modern battlefieli is comprised of many class2s and 
quantities of seismic signals. Using th2se signals for 
target identification and acquisition is the rurpose of 

re) 


| 
i 


batclefield seismic sensors. AS presentei in the last 
Chapter, matched filters can be used for target iderntifica- 


tion. In this chapter, it will be shown chat matched filter 


=] 


information may also be used to obtain target bearing. 

Multiple target acquisition raquires she concurrent 
separation of tha target classes and the computation of 
direction for each of the target classas found. The output 
of the matched filter is a spixe at t~, the cine of peak 
signal detection. If matched filtering*eis performed for 
each of the sensors in the ring , the value of «heir respec- 
tive t,, for each filter output, woild be different. Since 
+he seismic waves impinge upon each sensor at different 
times, depenient upen the target diracticn, the values for 
noe for each sensor may be expected to be directly ralated 
to the arrival angle of the seismic wave. The ¢, spike of 
the matched filter output may be thought of as a_ signal 
compression for hoth the continuous (tank, truck etc.) and 
the time limited (shell blast, artillery recoil etc.) 
seismic signals (Ref. 6]. 

Time domain aethcds (Ref. 1] use the positional differ- 
ences of known wave points to geomstrically estimate direc- 


tion. The tiaes associated with a sensor ring*s peak 
amplitude responses may be explained by way of an illustra- 
tive example of a shell blast. A rough direction to the 


origin of this shell blast can be computed using the 
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relative tine differences associated with ‘he peak 
amplitudes of all of the ring's sensors (Ref. 9]. 


oy 


The enhanced time positional information, which is 
by-product of the matched Filtering, san be used te perfora 
just such a time domain approach. The motivation for the 
method is that, unlike other cine and frequency domain 
acthods, which are very susceptible <co noise ccrruptio 


re | 


inaccuracies, matched filtering pulls the signal out of the 
noise and optimally detects she signal at time t.. Two tims 
dcmain methods are evaluated for finding arrival angles of 
the seismic signal. The first being the time domain chase 
difference (TOPD) snethod. The s2cond being a leas* nean 
squares polynomial (LMSP) curve fitting approach. 


B. MULTIPLE TARGET PILTERING ALGORITHS (MULTI) 


The routine MULTI calls the match2d filter routine for 
each sensor's amplitude signal. Returned are the target 
classes found with their relative peak filter response posi- 
tioas. Figure 5.1 illustrates a two target case. The 
matched filter response and the relative time postion for 
the two classes of targets can be seen. The figure shows 
that a shell blast target is present. The relative time for 
this target class is 5800. The ralative time returned for 
the simulated wheeled vehicle target is 3000. These tines 
er positions are relative Since each sensor's filter 
rasponse peaks are offset in time with respect to the peaks 
of the other sensors. This allows for simultaneous direc- 
tion finding for each class of target. This algoritha 
allows selection of either a time domain phase difference or 
or a least means square polynosial algoritha for finding the 
Girection to the targets. The time domain phase difference 
algorithm is derived first [Ref. 1]. The least mean squares 
polynemial direction algorithm then follows and is believed 
+o be an orignial application to this field [{RBef. 10}. 
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pendence. 
2 (I - INA 


Multiple Target Direction Phase Difference Aigorssin 


Tj 7° the arrival time of the saisaic signai at 
a 


the I-th sensor for the J-th target class 


w 
tad 
( 


a - the angle of the I-th sensor to th 


Dj 7 the distance from the origin to where the 
wave front of the J-th target passes the I-th 


sensor 


To Ot the time when the J-ta target class passes 


the orgin 
Ki ,% - the position of the I-th sensor 


B,; - the arrival angla of the seismic wave for 


2 


the J-th target class 
Vo the seismic wav2 veiocity 
R ~ the radius of the sensor ting 


I - ‘the sensor number where I has integer values 


form one to nine 
J - the J-th target class 


N - the number of sansors in the ring 


Figure 5.2 illustrates these parameters and “heir interde- 


The derivation of the algorithm follows: 


where zero degrees is set parzallel to the x - axis 
» ad Rcos®@; 


= Rsind; 


= Rcoos(i - B ) 


| 
| 


Figure 5.2 Circular Sensor Array Geoaetry 


| or equivalently D.. = (X.)cosB_ + (Y )cosB. 
| ji i J i a 
Te, = ayy 2, 
i=l 
T54 = Te, 7 T5i 
ye, Dery: 
v= -)cosB. + (Y¥. )sinB. Te. -T.. 5.1 
((%) 5 (% ) 7k 5 3a? ( ) 


Since the wave velocity can be assumed to be 
constant when passing all sensors, thea for any sensor I and 
K where I # K, equation 5.1 leads to 


((X;) cosB, + (%, )sinB;)/ (Te; - 


T 
((X,) cosB. + CF) sipB, VATE. - Be 


where I # K 


Now solving for the arrivai angie 2. 
J 


= - - (Ts = : } 
Bik arctan ( (is Te ) x. ( Ti )X )7 
((To =~ T..)Y¥ - (7e - T+)Y~)) 
3 ji ok J 


or equivalently; 


Bs ix = arctaa( ((To; -T.,)cosd - (Tc, -T,,) cosh )/ 


((fo, -T, y sind (Te -T) S228. 2) 


Where thé values of T.; and T., are returce values frtm tne 


matched filter routine. 


Now; 


: s 2 
B; C14 (N) Ere 
L=1lk=1l 
Whare 8B; is the direction in cadians to <+he J-*h class 
n 


target. For the multiple direction routine as i 
'j' has values from cne tc five. 


2. Jeast Mean Square Polynomial Direction Finding 


The least mean square direction finding aigeoritha 
was developed in response to probiems encountered with the 
phase differance direction finding algozitiua. This xew 
method is base on a ieast mean squares polynomial curve fit 
of «he sensor data. This approach was selected since «he 
least mean sguares polynomial provides for best fit oz 4a 
maxioum likelihood curve fit for noisy data. 

The laast mean squares direction finding algoritha, 
as with the phase difference algorithm, assumes the seismic 
wave to be planar. Pigure 5.2 illustrates the parameters 
for this model. Once the assumption of a plannar seisnic 
wave is made, the expected relation between the arrival 


angle, relative delay times 
circular array, can be nada, 
relations fcr a nine sensor ci 
. Wave arriving at zero degrees. 
delay «ines have been scaled to 
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Pigure 5.3 Relative Delay Times in a Nine Sensor Ring 


pie Seals Sn 


Since the sensors in a 
inerements cf forty degrees relative to the 
lation can be seen bétween the rei 
sensor and the arrival angle of cn2 s@ismic wav 

a 


can now be made <¢to illustrate th2 ralationship of san 
angles versus relative time delays. Pigura 5.3 is a plot 
for a nine sensor ring with a plansnar seismic wave arrivin 


at zero degrees. Figur)? 5.4 shows that for 2 wave at an 
arrival angl2 of 160 degrees, cane delay ‘ine at sensor 


ct 
+ 
rh) 
th 
Ve 
oF 
it 


number five will be zero. I* can be seen tha 
of these ideal data points with a l2ast mean sq 
mial will produce an equation for a curve wh 
value is also at the arrival angle 9f the seisn 
Minimum degree of ‘the pelynomial ts fit this i 
four. This results from noting that ¢2 
has three curve inflections. For expsrinent a 
minimum curve point corresponds to the predicted arri 
angle. : 

Polynomials cf degrees higher thar four aay be 
expected to enhance arrival angle. errors since tha pelyno- 
mial would distort te fit noisy data. Least neans squares 
polynomials of degree two and chres however, may be useful 
in reducing curve sensitivity to one or two aalfunctionin 


sensors or excessively noisy data. 
3. Least Sean Squares Poiynomial Algoriths Derivetion 


Let: 
N - number of data point pairs 


Y, - the observed or axpetimatal data position 


values 


Xi- independent degree values with a range of 0 to 
360 degrees 
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Figure 5.4 Relative Time Delay versus Sensor Angle 


E> dependent predicted dalay time values found by 
the least mean squares polynomial 


A, -B; - coefficient values for the system of simul- 
taneous equations 


E,~ error between the experimental data values and 
predicted values of time dialays 


S ~ sum of the square errors betweer each data and 
predicted points 


The derivation of the least mean squares polynomial follows: 


S$ = £2 ¢ E2 + £2 #..+ 52 
1 2 3 N 
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n 


P= A+ ALK + ASX? + AKS # oo. AK! (5.2) 

Be ¥ 7 Py (52) 

s = Se (5.4) 
i=l 


Combining equations 5.2 and 5.3 yields 


n 
Z£.= b ies Ay,~ & Xin AX? ~~ eee” 4a, & 


eae U L 


where n is the degree of the polynomiai such that N> net 
ard 1 <i<N 

Equation 5.4 now becomes, aftsr substituting for £, equation 
5.5 


S = vy, - ee See Ve mat ia A xi2 (5.5) 
i=l 
T> find the minimum cf the sum of tha sdgduares express¢i by 
equation 5.5, the partial derivatives of S with respect t9 
all of the coefficients are <axken. At the mininum, these 


partial derivatives all vanish. 
N 
YS/A_z 0 = 22 (Wir Age A Myr eee AKG) (OD 
L1= 


BS/9A = 0 = D2 Rom By Kyron ee A, Xi) (CX) 


8s/da,2 0 = S2(%p- Aom A Xyreee- AEG) (-FE) 


Dividing by two and rearranging gives fc + 1 normal simulta- 
neous equations. Expressed in matrix notation these squa- 


tions becomss 


yn dx, x2... Fx; {fa 
2K xe Dks we. Dx, [fa 
dx x3 Dxe es De. iia 


rr Yyntl Yana. Dxen tia Deer | 

The coefficient matrix in the abcve system of equations can 
be solved for. The mininum value for Pp. can then be found. 
The X corresponding te this minimum value is declared *o be 
the predicted arrival angle of the seismic wave. (Ref. 10] 


4. Adaptive Target Direction Finding 


An adaptive method is used to inprove the directicns 
found for all single target cases. To enhance the resolu- 
tion accuracy of the peak filter output position, the single 
target seismic data is copied into th2 filter data section 
for its class of target. Subsequently, when the matched 
filter routine is called for each sensor, filtering is 
performed with sensor data from that time period. 


S- Software features of the Sultiple Zarget Direction 


As first addressed in tha previous chapter, the 
multiple direction routine allows for the selection of the 
number of data points to be used from the 1024 size buffer 
of sensor dati. This was included to investigate the algo- 
rithe's performance for various data/filter window sizes and 


to allow an option for reduced <CPU time utilzation for 
interactive program runs. Tabular output from this routine 
includes the event number, the time period for the data, 
identification and directions for up to five target classes, 
and up to four sisulated target spacifications. Also, 
notice that the routine adapts t)5 the number of sensors 


specified for the ring. This was nacessary since the exper- 
imental data included sensor rings of six and nine ver-ical 


Sensors. 


C. SULTIPLE TARGETS OF THE SAME TARGET CLASS 


A limitation on the system, is presented up to this 
point, is its inablity t9 engage multiple targets of the 
same class. Discrete targets may appear as separate enti- 
ties since the seismic signals ganerated by <he two same 
class targets area not likely to be incident at the same 
time. This situation is greatly complicated for continucus 
tine signals, such as tracked or whaelead vehicles. For such 
continuous time targets, mutual distortion would be the 
likely result. 

An algerithm is needed to identify these multiple peaks 
without erroneously declaring elemsnts in the same peak as 
targets. By using cnly the positive half of the matched 
filter output and performing data smoothing on the remaining 
points, a curve with its number of peaks equalling <¢he 
number of targets present could be generated. Numerical 
methods for curve fitting or intarpolation are available 
(Ref. W1.J. Polynomial curve fitting, Sterling's method and 
variations of Newton's method are oaly a few of the possible 
approaches applicable for equally spaced data. By differen- 
tiating the smoothed data, the paaks and vwalleys of the 
filter output can be found. The height of the curve corre- 
sponding to the points where the dsarivative is zero can now 
be compared to the selectad matched filter threshold. This 
excludes the valley foints and leaves the points remaining 
which correspond to the relative tines of the same class 
targets. 

This method is ccnstrained by resolution of close prorx- 


imity, time limited targets and aaar phase synchronization 


ef continuous targets. The variables 
through experimentation may be the degree 
of the curve data andthe exclusion of 


associated with the same target's data. 


to «6be)| 6 6optinized 
of the smoothing 
erronecus valleys 
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VI. ANALYSES OF SEIZSHIC DATA 


Analysis of the sinulated and 2xperimental seismic data 
will be conducted as detailed in Tabie II. Table ITI lists 
the matched filter contents for th2 simulated data. Table 
Iv is the test plan for the experimental data. Table V 
lists the matched filter signals used for the experimental 
data analysis. 

Table VI summarizes the results of the simulted and 
experimental data runs for direction fiading. The window 
size, used for all multiple cscarget direction finding 
results, was 300. Table VII lists the events in which 
targets were nissed or incorrectly identified. 

The time domain phase differenc2 directions found, are 
not presented for the reasons noted earlier. Errors cf up 
to aighty degrees were not uncommon with this method. 

Bach event run will be accompanied by the follcwing 
graphic output: 

1. Least 4ean Square Initial Direction 

2. Matched Pilter Response 

3. Amplitude Respcnse 

4. Amplitude Response of any Malfunctioning Sensor 

5. Frequency 2esrfonse 

6. Least Mean Squares Polynomial Curve Pitting (LMSP) 
Using Matched Filter Outputs 

7. j\Multiple Target Direction Summary Resulting from 
Least Mean Squares Curve Fitting 


a7 


TABLE II 
Test Plan for Simulated Data 


10 30090 40 
10 3000 129 
10 3000 240 


oO 

o 

- 
oOo © 0 
~_ 2 = = 


2 
\ 
f 

| | 
Event #sen f#Iagts Frequency Amplitude Jzzectson 

| 10 3090 9 | 
| \ 
| 
| 

| 

| | 
4 


TABLE IIIf 
Matched Filter for Simulated Targets 


Filter Preguency Amplitude Direction 
Tracked Veh 30 2930 0 
Wheeled Veh 10 2090 ¢] 
Shell blast data from event #333 4) 
Helicopter 15 2020 0 
Personnel 20 2090 0 


EE ea caret Ae oe a A Oe ee a eee 


eA AAD pete OS ag eee OE GO SO een SS tae eiRanty 


Eee etter one RE- EIEIO SES SE TSE men mseeerersnece asses 


_ SS SS TT SS TT SS aS a 


TABLE IV 
Test Plan for Experizantal Data 


| 

! 

| | 

{ Evens #sen #fats Dir Target Distance 
383 9 1 0 Shot SKM | | 
382 9 1 0 Shot 5kn { 
372 6 1 315 Helicopter 5 - 15K | 
375 6 1 0 Taak 5 - OKo { 

374 6 1 315 Helicopter 15 - SKM 

| 302 6 1 0 Mortar 1K ! 

« #314 6 1 315 Lvt 4 - SKM { 

354 6 5 0 105mm How 5K | 

225 175mm Gun 4KM ! 

315 LvVtT & = SKM { 

| 0) M-60 Tank 4& = SKM 


ee 


TABLE V 
Matched Filter for Experciaental Data 


| Target Event Used as Piltes 
Tracked Vehicle 375 

Blast;sRecoil 383 

Helicopter 372 


Personnel none (backround noise) 


{ 
| 
{ 
| 
| 
Wheeled Vehicle Rone (backround noiss) 
| 
( 
| 


 - 


ee =. ~q 
| 
| TABLE VI | 
Summary of Diraction Pinding Results | 
, doptial LMSP | 
Zver: #Sen #Igts Diszance Diz DLE Zezror Dit €>ECz | 
| 001 9 1 N/A 0 N/A ) 0 | , 
1 001 9 1 N/A 4QO N/A 40 0 
| 901 9 | N/A 120 N/A 120 0O 
001 9 1 N/A 240 N/A 240 O 
383 9 1 SKM 9 4(3) 1.1 28(4) 7.8 | 
{ 5 1.4 
| 382 9 1 5KM 0 -14(4 3.9 -6 1.67 
372 6 1 S - 15K 315 ~32 3.6 -59 3.9 | 
{ 375 6 1 5 - OKM QO 312 13.3 PAILED | ie 
| 374 6 1 15 - SKM 315 -59 4.0 FAILED i. 
! 302 6 1 1KM 0 4 1.1 6 1.7 
354 6 5 5KM 0 >3 -8 3 8 { 
4 KM 225 291 18.0 PAILED 
G4 =-5KM 315 FAILED 
4 - 5KM O PAILED 
*Note: Tha matched filter threshold was | 
| at_.9 for all single taryets and 6 “for all | 
multiple target events. Srackets indicate the 
us2 of other than a sa3cond degree polynomial. 
! { 
ene = See 
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pase wake numbered error sources apply to all events 
listed. 


eres - (2 
| 
| TABLE VII 
| Sissed or Incorrectly Identified Targets 
Event #Targets Target Nature of Errors | 
: { 

| 375 1 Tank 1. Matched filter was not i 

pased on a high S/N 

| sample signal j 

374 1 Helo 2. Small seismic signal | 

| anplitudes ; 

354 1 LVT. 3. Malfunctioning sensor(s) | 

| 1 Tank aude cts : : \ 
( 1 175am Sever ciipping distortion of| 
| Gun anput signal 
( 
{ 
t 
{ ! 
{ | 
4 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


The ability of the digital matched filter to dexect anal 
correctly identify actual discrete single targe: types was 
excellent. The adaptive anhancement of the matched 
schame was found to sharpen the matched filter responses. 
The complications of multiple targets and continuous cargets 
proved to be less successful. Lack of high signal to noise 
ratio sample signals for use as filters reduced the apileiy 


of the filters to match the signals. As would be 


Oo 


xpected, 
the Simulated target identification and direction finding 
Operations met with success for poth single and multipl2 
targets. 

The simulated data validatsd the usefulness of the least 
mean squares curve fitting methoi ‘for target direction 
finding. This algorithm was noted to be useful in both tha 
relative peak amplitude response method for recoily/blast 
targets and the matched filter peak position method for all 
targets. The highest accurracies were found using second 
degree polynoaials. This was due to the reduced noise 
sensitivity of lower degrae poiynomials. Cconvecsly, signif- 
icant errors were found in ‘the diractions determined by the 
phase difference algcorithe. Thesa errors were possibly due 
to round-off error sensitivity is the software/hardware 
inmplesentation. Experimental data could not be used to 
effectively crosscheck this finding since most of the exper- 
imental data targets were at zeros degrees and the phase 
difference routine seemed to seek zars degrees. 

“a significant resnlt was the accuracy of the bilast/ 
tecoil target directicn found using ouly the peak asplitude 
responses. Directions could be found using only the rela- 
tive amplitude peak positions for the array sensors and the 
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least mean squares curve fitting routine. This resul- indi- 


cated a possible counter-fire applicaticn using a vrseclv 
simplified systen. This finding is felt +o be sign 
ince artillery type targets are the highest prioricy taryex 
type. The discrete nature of the blast/racoil seisa 
signals wouli also allow for ready separation of sv 
large number of ccmbined hostile and friendly signal 
sources. This would allow for observerless adjustment of 
fire onto hostile targets. Artillery and mcrtar targets 
may, infact, ba the only <target types detectable at the 
Tanges specified for a long tange seismic system. 


It is recommeded that further study be made in to the 


u 
possible implementation of the least mean square curve 
fitting of the peak sensor amplitudes responses in a 
counter-fire system. Digital matched filters for seismic 
target identification ani direction finding may also orove 
effective after further experimentation with optimum sample 
Signal filters for the various targat «ypes. Additionally, 
the matched filtar response may possibly be enhanced by 
preprocessirg the seismic signals through adaptive noise 
cancellors. 

Acoustic vice seismic matchad filtaring with directicnal 
microphones may be useful in target identification. Once 


the target has been identified, a aatched filter/leas* mean 
Sguares based direction finding scheme could then be 
attampted. 


ARPENDIX Aa 


USERS MANUAL 


The software developed provided for interactive progran 
oparation. However, further inforcmatien must be or 
fer an initial system setup and correct program operation. 
To begin with, the seismic data aust be transferred from 
Magnatic tave to the IBM 3033's Mass Storaq2 System (MSS). 
The data must then be «ransfered ‘+9 the DISSPLA user's disk 
for analysis of the program. These data transfers may be 
accomplished using Jcb Control Language (JCL) procedures. 
The magnetic tape volume must first be scannad to deter- 
mine the storage format of its files. The JCL procedure 
TSCAN provides this information. A sample TSCAN job 
follows: 


//ILIVIETT JOB (3026,0304),*SMC-1677 JOHNSTON',CLASS=F 
// EXEC TSCAN,VOLIN=PARKI, DCBIN=* DEN=2', UNITIN=a"Z 400-48 
// EXEC TSCAN, VOLIN=PARK2, DCBIN="DEN=2* ,UNITIN='3400-4' 
// EXEC TSCAN,VOLIN=PARK3, DCBIN=" DEN=2', UNITIN=" 3400-4" 
7/ EXEC TSCAN, VOLIN=PARK4, DCBIN="DEN=2*, UNITIN=* 3400-48 
1 


Once the tape scan is coaplated, the tape files and 
comments can be transferred to the SSS. Prier to this 
transfer however, space inthe MSS must be made for these 
files. The procedure ITEFBR14 is used for this purpose. A 
saaple jeb follows: 


//TLITSTGA JOB (3026,0308) ,*S8C1677 JOBNSTON' ,CLASS=A 
//*SEAIN ORGAHPG VMI. 0 131P 

//- =—-EXEC PGMSLEFBR14 

7/001 DO ONITS3330V, ASVGP=PUBGB,DISP= (NEW, CAT LG) , 

4/ SPACE (CYL, (4,4,3) ) » DSHS8SS. S30 26.P302 
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/* 

// EXEC PGM=LEFSR14 

7/001 DD UNIT#=3330V, MSVGP=P UBYB,DISP= (NEW,CATLS), 
// SPACE= (CYL, (4,4,3) ), DSN=MSS.S3026.P3 14 
/* 

// SXEC PGM=IEPBRI4 

//0D1 DD UNIT=3330V, MSVGP=PUBY4B,DISP=(NEW,CATLS), 
4/ SPACE= (CYL, (4,4,3) ), DSN=MSS.530 26.7319 
f* 

// 


Each event is proceeded by a comment ‘file for that 
event. These comment files can be identified from the TSCAN 
output as a file containing only one record. The files 
containing nine or eighteen recoris are the sensor data 
files for the events. Files of nine records in length are 
evants using 21 circular array of nine sensors designated as 


a type A33 array. These sensors are all vertical aoticn- 
sensing geophones. The 2ighteen record files contain six 
senscr groups of three geophones. A circular érray is 


designated type A31 and a linear array is a type A32 array. 
The three geophones for each group sense either -cradial, 
transverse cr vertical motion. A JCL routine, using the 
procedure IEBGENER, transfers the event comment and sensor 
Gata. A sample IEBGENNER job follows: 


J/ILIVIGT7T = JOB (3026,0304), "JOHNSTON SMC1677" ,CLASS#P 
7/*4RIN ORG=NPGVM1.0131P 


//* 

//* CRE/CHS SUBMIT TEBGENER JCL 

s/* 

s/* COPY TAPE PILES TO 8SS.S$30 26.7319 
//* 


// BXBC PGHELEBGENES 
//SUSPRINT DD SYSOUT#A 
//STSUT1 DD UNIT#3 400-4, VIL =SER=PARK 2, DISP®(, PASS), 


119 


// LABEL=(1,NL,,1IN), 
// + DCB=(RECPM=FB, LRECL=64 , BLKSIZE=2 948, DEN=2, IPTCD=0) 
//SLSUT2 DD DISP=SHR,DSN=4 SS. $3025. COMMENTS (COM3 19) 
//SYSIN DD DUMMY 

f//* * &€ ee ee eK KKK SE RK KK KEK KK KF 

//COPY PROC FILE=,MEM= 

// EXEC PGM=IEBGENER | 
//STSPRINT DD SYSOUT=A 
//SYSOUTV DD UNIT=3400-4 ,VOL=SER=PARK2,DISP=(,PASS), 

1 LABEL=(&FILE,NL,,IN) ,DCB= (RECFA= PF, BLKSIZ2=2048, DEN=2) 

//SYSUT2 DD DISP=SHH,DSN=MSS.$3026.2319 (EMEM) 

//SXSIN DD DUMMY 
// *PEND 

/* 

Sf/* * 8 eR KE KEE KH KK KF HK RK KK 
//—«=EXEC COPY, PILE=342,MEM=SEN1 

//-BXEC COPY, FILE=343, MEN=SEN2 

//—«BXEC COPY, PILEZ=344, MEM =SEN3 

//—«OSXEC COPY, PILE=345, MEM=SENG 

//- «EXEC COPY, PILZ=346, MEM=SENS 

//—~xEX8C COPY, FILE=347, MEN=SEN6 

// —-«BXEC COPY, PILE=348, MEM=SEN7 

// _—*-BXEC COPY, PILE=349, MEM=SENS 

//—«xBXEC COPY, PILE=350, MEN=SEN9 

//—s«xBXBC COPY, FILE=351, MEM=SEN10 

// EXEC COPY, PILE=352,MEM=SENT1 

//-—«xBXEC COPY, PILE=353, MEM=SEN12 

//—*xBXEC COPY, PILE=354, MEM=SEN13 

// ‘EXEC COPY, PILE= 355, MEM=SEN14 

// —*xBXEC COPY, FILE =356, MEM =SEN1S 

//-—sOBXEBC COPY, PILE= 357, SEN=SENI6 

//—« EXEC COPY, PILE=358, MEN=SEN17 

//—«BXBC COPY, FILE2359, MEN=SEN18 

//_—*xBYEC COPY, PILE2360, SES#SEN19 
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// 


Transfer of the ¢ ::>22 avens data from SS <= “h2 
DISSPLA user's disk may .cw °°. performad. A patch fortran 
job with the appropr- LLEDEFPs +9 déenote the varicus 


geophone's data is submitted. The RSC 
te send the output of this routine +9 


saaple fortran job for nin2 sensors follows: 


7/314983026 JOB (3026,0304),' JOHNSTON’ ,CLASS=a 
//* MAIN ORG=NPGVN1.0090P 

// EXEC FOFTXCG,REGION.GO=1024K 

7/FORT.SYSIN OD ® 


Cc 
LOGICAL*®1 INFO1(8) ,IN FO2 (8) , =INFO3 (3) - INFO4 (8), INFOS (8) 
LOGICAL*1 IN?0O6 (8) ,INFOT (8) ,IN?F03 (3) , INPOI(8) 
INTEGER*2 DATA1(1020) ,DATA2Z (1229) ,DATAI (1920) , DATAY (1020) 
INTEGER®*2 DATAS (1020) , DATA6 (1920) , DATAT (1020) , DATAB (1020) 
INTEGER*2 DATA9 (1020) ,DATB1 (4990) 

Cc 


DO 30 J=1,4 


READ(1, 100) INFO1, DATAI 
READ(2, 100) INFO2,DATA2 
READ (3, 100) INFO3,DATA3 
READ(4, 10.0) INFOU,DATAY 
READ(8, 109) INFOS, DATAS 
READ(9, 100) INPO6, DATAG 
READ(10,100) INFO7, DAT AT 
READ(11,100) INFOS, DATAS 
READ(12,100) INFOS, DATAS 
100 PORHAT (8A 1, 10.2 (1082) ) 
DO 10 fT = 10, 1020,10 
WRITE (6, 701) DATAV(I ~ 9) ,DATAT(I - 8) ,DATAU(I - 7), 
1IDATAT (I - 6),DATAT(I - 5),DATAI(I - 4),DAPTAI(I - 3), 
2DaTA? (I - 2),DATAI(I - 1),DATA1 (I) 
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10 


20 


4Q 


50 


60 


70 


80 


o> =< wraene stays eee lis ce 5 


CONTINUE 

po 20 I = 10,1020,19 

WRITE (6, 101) DATA2(I ~ 9) ,-DATAZ(I ~ 8),DATAZ(T - 7), 
ADATAZ(I - 6) ,DATAZ(I - 5),DATAZ(I ~ 4),DATA2(I - 3)¢ 
2DATA2(I - 2 ,DATAZ(I - 1),DATA2 (I) 

CONTINUE 

po 40 I = 10,1020,10 

WRITE (6,101) DATAS(I - 9),DATAI(I - 8),DATA3(I - Te 
IDATA3Z(I - 6) -DATAZ(Z - 5),DATA3(I - 4),9ATA3(I - Fe 
QDATA3Z(I ~ 2),DATA3(I - 1), DATA3 (T) 

CONTINUE 

po 50 I = 10,1020,10 

WRITE (6,101) DATAG(I - 9) ,DATAS(I - 3),DATAY(I - 7), 
WDATAG (I - 6) pDATAY(I - 5),DATAS(I - 4),DATAU(T ~ 3) 
2DATAG(I - 2),DATAS(I - 1),DATAS (I) 

CONTINUE 

pO 60 I = 10,1020,10 

WRITE (6,101) DATAS (I - 9) ,DATAS(Z - 8),DATAS(I ~ 7), 
WATAS (I - 6) ,DATAS(I - 5),DATAS(I - 4),/DATAS(I ~ 3)» 
QDATAS(I - 2) ,DATAS(I - 1),DATAS (I) 

CONTI NUZ 

pO 70 I = 10,1020,10 

WRITE (6, 101) DATAG (I - 9) ,DATAS(Z - 8),DATAG(I ~ 7) 
WDATAG (I - 6) ,DATAG (I - 5),DATAG (I - 4), DATAG(I ~ 3) 
2DATA6 (I ~ 2),DATAG(I ~ 1),DATA6 (I) 

CONTINUE 

po 80 I = 10,1020,10 

WRITE (6,101) DATA7(I ~ 9) -DATA7(I - 8),DATA7(I - 7), 
IDATA7 (I - 6),DATAT(I - 5) ,DATA7(I ~ 4),DATAT(I ~ 3), 
2DATA7(I - 2,DATA7T(I - 1),DaTA7 (I) 

CONTI NUE 

po 90 I = 10,1020,10 

WRITE (6,101) DATAG(I = 9) ,DATAS(I - 8),DATAS(I ~ 7), 
(DATAG (ZL - 6) ,DATAS(I - 5),DATAS(I - 4), DATAB(I - 3)- 
QDATAG(I - 2),DATAS(I - 1),DaTA8 (I) 
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Fe a mgr oe 


90 CONTINUE 
DO 91 I = 10,1020,10 
WRITE (6,101) DATAI(I - 9) ,DATAI(I - 8),DATAQ(I - 7), 
I1DATA9 (I - 6) ,DATAS(I - 5),DATAO(I - 4),DATAQ(I - 3), 
2DATA9(I - 2),DATAQ(I - 1),DATA9 (I) 


91 CONTINUE 
101 FORMAT(I6,16,16,16,16,16,16,16,16,16) 
30 CONTINUE 
STOP 
END 
/* 


//30.FTOIFO01 DD DISP=SHR, DSN=%SS.S30 26.2383 (SEN1) 
//GO.FTO2F001 DD DISP=SHR,D SN=MSS.5S36 26. P383(SEN2) 
7/30. FT03F001 DD DIS P=SHR, DSN=4SS.530 26. P383 (SEN3) 
//GO.FTOWFO01 DD OILS P=SHR,DSN=MSS.S3026.7383(SEN4) 
//GO.FTO8FO0T DD DISP=SHR,DSN=MNSS.S30 26. F383 (SENS) 
//GO.FTO9PO0T DD DIS P=SHR, DSN=NSS.S30 26. P383 (SEN6) 
//GO.FTIOFO91 DD DIS P=SHR,D SN=MSS.S30 26. P383(SEN7) 
//30.FTUIFOO1 DD DIS P=SHR,DSN=MSS.530 26. P3833 (SEN8) 
7/GO.FT12F001 OD DIS P=SHR, DSN=NSS.S30 26.2383 (SEN9) 
7/GO.SYSIN DD * 

/* 

// 


A sample job for a six sensor group array follows: 


//3L383026 JOB (3026,0304), ‘JOHNSTON',CLASS=A 

//*4AIN ORG=NPGVM1.0C90P 

// ®XEC PORTXCG,REG ION. GO= i024K 

//PORT.SYSIN DD * 

c 
LOGICAL®1 INPO1(8) ,INPO2 (8) ,INPO3 (8) , INFOU (8} , INPOS (8) 
LOGICAL*®1 INPO6 (8) , IW PO? (8) , INFOS (8) , INFO9 (8) 
INTEGER*2 DATA1(1020) ,DATA2(1220) , DATA3 (1020) , DATA (1020) 
INTEGER®2 DATAS (1020) -DATAG (1220) ,DATA? (1020) , DATAS (10 20) 
IBTEGER*2 DATA9 (1020) , DATB1 (4996) 
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DO 30 Jz1,4 


READ(1, 100) INFO1, DATA 
READ(2, 100) INFO2,DATA2 
READ(3, 100) INFO3, DATA3 
READ (4, 100) INFOU,DATAY 
READ(8, 100) INFOS, DATAS 
READ(9, 100) INFO6,DATA6 
100 FORMAT (8A 1,102 (1042) ) 
po 10 r = 10,1020,10 
WRITE (6,101) DATAI(T - 9) ,DATAI(I - By, DATAI(I - Ve 
IWATAT(I - 6) ,DATAIT(I ~- 5), DATAI(I - ¥),DATAT(IT ~ 3), 
QDATAI(I - 2),DATAU(I ~- 1),DATAT (TL) 
10 CONTINUE 
po 20 r = 10,1020, 10 
WRITE (6, 101) DATAZ (I -~ 9) ,DATAZ(I - 8),DATA2(I - TV)e 
IDATAZ(I - 6) ,DATAZ(I - 5),DATAZ(Z - 4), DATA2(T ~ ee 
OpATA2(I - 2),DATAZ(I ~ 1),DATA2 (I) 
20 CONTINUE 
po 40 r = 10,1020,10 
WRITE (6,101) DATAI(I ~ 9) ,DATA3(I - 8) -DATAZ(T ~ 1), 
wWatazs(I - 6) ,DATA3(I - 5), DATA3(I - 4), DATAI(I ~ 3). 
2DATA3(I - 2),DATA3Z(I - 1),DATA3 (1) 
40 CONTINUE 
po SO I = 10,1020,10 
WRITE (6, 101) DATAG (I - 9) ,DATAY(I ~ 8) -DATAS(E ~ Ve 
WATAG (IL - 6) ,DATAY (I - 5),DATA4(I ~ 4), DATAS(T ~ 3), 
QDATAU (I - 2),DATAM(I ~- 1),DATA4 (I) 
50 CONTINUE 
po 60 © = 10,1020, 10 
BRITE (6, 101) DATAS (I ~ 9) -DATAS(I - 8),DATA5S(T ~ 7), 
IDATAS (I - 6) -DATAS(I - 5) ,DATAS(Z ~ 4),DATAS(T ~ 3) 
2pATAS (I - 2 ,DATAS(I - 1)-DaTAS (I) 
60 CONTI SUE 


san 
ral ena casi 
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DO 70 I = 10,1020, 10 
WRITE (6,101) DATAG (I - 9) ,DATAS (I ~ 8) ,DATA6 (I - 7), 
IDATA6 (I ~- 6) ,DATA6(I - 5),DATAS (I - 4),0ATA6 (I - 3), 
2DATA6 (I - 2),DATAG(I - 1), DATA6 (TI) 

70 CONTINUE 


101 FORMAT (I6,16,16,16,16,16,16,16,16,16) 
30 CONTINUE 

STOP 

END 
/* 


//GO. FTO 1F001 DD DIS P=SHR,D SN=MSS.S30 26. P350 (SEN3) 

//GO.FTOZFO0T DD DIS P=SHR, DSN=NSS.S3026.7350(SEN6) 

//GO.FTO3FO001 DD DISP=SHR,DSN=MSS.530 26. P350(SEN9) 

//GO.FTOSFOOT DD DISP=SHR, DSN=MSS.S30 26. P350 (SEN12) 
//30.FTO8FO0T DD DIS F=SHR, DSN=MSS.S3026.P350(SEN15) 
//G0.FTO9F001 DD OISP=SHR,DSN=MSS.53026. P350 (SEN18) 
7/GO.SYSIN DD * 

/* 

/f 


Two files will be returned to user's reader. fhe first file 
is the listing and diagnostics file and should be purged. 
The second file should be named SEN DATA. SEN DATA must now 
be edited. Delete the first seven lines of the file and 
issue the coamand LREC 80 to set th3 proper file record 
length. 

The interactive program can be run with the completes | 
collection of files listed in the exec MATCH. The MATCH 
ExzC follows: 


STRACE OFF 

FORTGI SFILTER 

GLOBAL TXTLIB PORTHOD2 SOD2EEH IMSLSPF VYONIMSL 
PILEDEF 10 TERMINAL 

FILEDEF 05 DISK SEN DATA (PERM) 

FILEDEY OF DISK COM DATA (PERS) 
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FI 4 DISK FILTER DATA (RECFM VS PERM 
FILEDEF 18 DISK DISSPLA METAFILE T4 (RECFPM VBS LRECL 
19065 BLOCK 19669 


EXEC DISSPLA MFILTER 


Entry parameters, such as the number cf sensors in *he 


array and the sampling rate, can be found in the NOSC data 
log for the event under study. Ali other interactive entries 
ara user salected options or are self explanatory. 
(Ref. 12) 


APPENDIX 8B 
SAMPLE INTERACTIVE PROSRAM SESSION 


match 
G1 COMPILER ENTERED 
SOURCE ANALYZED 
PROGRAM NAME = MAIN 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = ANGLE 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = TIMOUT 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = FREQOT 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = MATCH 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = MYDATA 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM NAME = MAXMIN 
* NO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAM WAME = SPCTRM 
: * WO DIAGNOSTICS GENERATED 
SOURCE ANALYZED 
PROGRAMS MAME = MULTI 
* NO DIAGROSTICS GENERATED 
SOURCE ANALYZED 


tk gerne ener 


PROGRAMS NAME = MLTPLT 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAM NAME = RMS 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAM NAME = AVG 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAM NAME = SIMULT 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAS NAME = LMS 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAM NAME = PLT 

* NO DIAGNOSTICS GENERATED 

SOURCE ANALYZED 

PROGRAM NAME = SOLV 

* NO DIAGNOSTICS GENERATED 
*STATISTICS* NO DIAGNOSTICS THIS STEP 

DISK *T* NOT ACCESSED. 

B (126) R/O 

C (127) R/O 

E (128) R/O 


ee. Your Portran program is now being loaded ... 


eee execution will soon follow ... 
EXZCOTION BEGINS... 


ENTER FVENT RECORDING NUNBER-I3- 
383 


EVENT NOMBER 383 


383. A33. 10 SEPTG1 SKE EOD SHOT. SE? B. BO DELAY. C141 CH FINAL 


ENTER 


ENTER 
1000. 


ENTER 
0400 


END OF TAPE, “ 


SAMPLE RATE IN HERTZ-REAL—~ 


LOW LOOK ANGLE IN DEGREES-I3- 


HIGH LOOK ANGLE IN DEGREES-I3- 


MATCH FILTER THRESHOLD. RANGE OF 0. TO 1.9 


PLOT SCALING FACTOR~ REAL- 


NOMBER OF SENSORS IN RING 6 OR 9 ONLY-I1- 


SENSOR NUMBER FOR DISPLAY-I1- 


NOISE THRESHOLD LEVEL -REAL- 


DATA WINDOW SIZE FOR DIRECTION FINDING-I4- 


FOR COMPRS OUTPUT ENTER -1-. FOR TEK618 ENTER -2- 


1 


TO CREATE SIMULATED TARGETS ENTER ~1- ELSE -2- 


1 


PUTER 


EUTER 
0 


ENTER 
0 


THE POUR SIAULATION FREQUENCIES-RZal- 


PREQUENCY 


PREQUENCY 


ENTER FREQUENCY 


ENTER FREQUENCY 


ENTER AMPLITUDES FOR EACH FREQUENCY -82AL~- 


ENTER AMPLITUDE 


ENTER AMPLITUDE 


ENTER ASPLITUDE 
0 


ENTER AMPLITUDE 
4000. 


ENTER TARGET ANGLE FOR EACH FREQUENCY-I3- 
SIX SENSORS ALLOWABLE ANGLES; 0,60,120, 180, 240,300 
POR NINE SENSORS;0,40,80,120,160,209,240, 280,320 


ENTER ANGLE 
0 


ENTER ANGLE 
0 


ENTER ANGLE 
0 


EATER ANGLE 
120 


TO SODIFY AMPLITUDE OF SIGNAL ABOVE NOISE THRESHOLD 


BNTER - 1 -,ELSE EWTER - 2 - 


130 


Oe En ae 


2 
ENTER DEGREE OF POLYNOMIAL DESIRED ~I1- 


FOR SIX SENSORS ENTER 2 - 4, FOR NINE ENTER 2 - 7 

4 

>> USING A PRE~ALLOCATED DATASET FOR UNIT FT17FO01. 
>> USING A PRE-ALLOCATED DATASET FIR UNIT FTISFOO1. 


| 

CATALOG AS TARG2I? IP YES TYP2 - 1 -, ELSE - 2 ~ | 
2 | 
4 


FOR MATCH FILTER DIRSC FINDING SNTER ~ 1 ~- SLSE -2- 
1 


FOR PHASE DELAY METHOD ENTER -1,FOR LAS ENTER -2 
2 
ENTER DEGREE OF POLYNOMIAL DESIRED -I1- 


FOR SIX SENSORS ENTER 2 - 4, FOR NINE ENTER 2- 7 
4 


TO VIEW OTHER SENSORS ENTER - 1 - 
1 


ENTER - 2 = IO CONTINUE TO NEXT TIME PRAMSE 
2 
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